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ABSTRACT 

A small neutrino Majorana mass can arise in the Standard Model as an 
effective dimension 5 operator. We calculate the renormalization of this 
operator in the minimal Standard Model and in its two-Higgs-doublet and 
super symmetric extensions. Renormalization from the scale of lepton number 
violation (e.g., the Planck scale or a GUT scale) to the weak scale decreases 
the strength of this operator by an order of magnitude or more if the top 
quark and Higgs boson masses are large. Neutrino mixing angles also run 
with momentum. We show instances where small mixing at a high scale 
becomes large at the weak scale, and vice versa. 
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1. Introduction. There are several indications that physical processes do 
not respect the global baryon number (B) and lepton number (L) symmetries. 
For example, in the Standard Model (SM), nonperturbative effects due to 
instantons break the B + L symmetry |IJ . The observed matter-antimatter 
asymmetry of the universe is an explicit evidence for B violation. In addition, 
simple arguments suggest that quantum effects of gravity may violate global 
symmetries 0. It thus appears that B and L symmetries are not fundamental 
symmetries of the full theory, but are "accidents" of the low energy structure 
of the theory — as is the case in Grand Unified Theories (GUT). If lepton 
number is indeed not conserved, a neutrino Majorana mass will be generated 
at some level. 

Let Mx be the scale above which lepton number is broken. This could be 
the Planck scale (if the lepton number violation is induced only by gravity) or 
a GUT scale. Assuming that the only light fields below Mx are the ones given 
in the SM,[] then the lowest dimension operator which generates a Majorana 
mass for the left-handed neutrinos, u Ll is unique JJ: 

C w = ^K mn l c ™ i ll j <& k <&ie ik e j i + h.c. (1) 

where II and $ are the left-handed lepton and Higgs boson doublets, respec- 
tively, 
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1 We do not consider here the possibility of light "exotic" particles such as triplet Higgs 
bosons, although later we shall consider the two-Higgs-doublet and the minimal super- 
symmetric extensions. 
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K mn is symmetric under interchange of m and n, the generation indices; 
i,j,k,l are SU(2) indices. When the Higgs scalar develops a vacuum ex- 
pectation value, ($°) = v/\/2, the SU(2) gauge symmetry is spontaneously 
broken and the neutrino Majorana mass matrix ensues: 

M v = ~kv 2 . (3) 

The operator in Eq. (|I|) has dimension 5, so k is of order 1/Mj. Conse- 
quently the magnitude of the neutrino masses is suppressed with respect to 
the charged fermion masses by the factor v/Mx- 

Models for neutrino masses predict the coefficients, n m n, in terms of other 
parameters. For example, if z/r's exist, they should have a large Majorana 
mass, Mr, from lepton number breaking and hence there will be a contribu- 
tion to k of 

K^2YjM^Y v , (4) 

which is the popular see-saw mechanism [||. Here Y v is the Yukawa coupling 
matrix between the l^s and i^rs. This prediction is not unique. Models with 
more than one heavy scale (such as the GUT and Planck scales), radiative 
models || || , or models with exotic heavy particles can give other predictions. 
The values of k depend on the short distance dynamics responsible for lepton 
number breaking, while the operator structure in Eq. (|l|) depends only on 
the low energy contents of the model. 

Radiative corrections renormalize the neutrino mass operator, as they 
do all the terms in the Lagrangian. Studies of the running of the SM cou- 
plings have led to the important observation that gauge coupling constants 
may unify at a large scale. The "running" of the neutrino mass operator 
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is also important for several reasons: 1) to correctly relate the scale of lep- 
ton number violation to other physical scales, e.g., proton decay, coupling 
constant unification, quantum gravity, etcQ and 2) to correctly relate the 
mixing angles predicted by a model at high scales to the values at low scales 
where they are measured. However, despite the fact that the renormaliza- 
tion of the dimension 6 operators which contribute to proton decay have been 
studied to two loops M, the running of the dimension 5 term which yields 
a neutrino mass has received little attention. In this note we compute the 
renormalization of the neutrino mass operator above the weak scale.f] We 
first consider the standard model with one Higgs doublet. Then we shall 
examine its two-Higgs-doublet and supersymmetric extensions. 

2. One Higgs Doublet. We choose to compute the renormalization of the 
neutrino mass operator, Eq. (|l]), in Landau gauge. The topologically distinct 
one-loop diagrams which give nonzero contributions are shown in Fig. 1. 
Beside the gauge interactions, there are the scalar quartic interactions, 

C H = ~(&$f, (5) 

and the Yukawa interactions, 

C Y = Ql^ur + Q L $Y}d R + l L $Yle R + h.c. (6) 

2 Solar neutrino experiments are sensitive to neutrino masses as small as 10~ 6 eV (see, 

e.g., [jjj). Such a small mass corresponds to lepton number breaking at a scale of order 

the Planck scale, k" 1 ~ ~ 10 19 GeV. 

3 The neutrino mass does not run significantly below the weak scale because there the 

neutrino interactions with other light particles are suppressed by the factor q 2 /M^. 
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Here Ql is the left-handed quark doublet, 



$ = zr 2 $* and ur, cIr, and cr are the right-handed quark and charged lepton 
fields. Yj is the Yukawa coupling matrix of fermion species / = u, d, e. 
Generation indices have been suppressed. With the above parametrization 
the physical Higgs boson mass is m 2 H = Xv 2 , where v = 246 GeV. 
The evolution equation for k is found to be 

Hk 1 
16VT 2 -— = [-3gl + 2X + 2S}k - -[n(Y^Y e ) + (tfY e ) T K ], (8) 
am/j 2 

where \i is the renormalization scale, g-i is the SU(2) gauge coupling constant 
and 

S = Tr[3Y^Y u + 3YjY d + Y*Y e ] (9) 

is the contribution associated with Fig. lc. 

The presence of two Higgs fields in the neutrino mass operator leads to 
qualitative and quantitative differences between the running of k and the 
running of the charged fermion Yukawa couplings 0. Firstly, A enters Eq. 
(H) at leading order, so the evolution is sensitive to large Higgs boson masses. 
Moreover, the evolution is twice as sensitive to a large top quark mass. Both 
of the above terms enter Eq. (|8|) with the same sign, so k can run faster than 
the charged fermion Yukawa couplings do. 

Solving Eq. (§), Fig. 2 plots how the magnitude of k evolves through 
the assumed desert above the weak scale. We numerically evolve all relevant 
coupling constants simultaneously over the entire range of fi between Mz 
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(the Z boson mass) and Mx, taken here to be near the Planck scale. (The 
renormalization group equations for the SM couplings can be found in, e.g., 



[y, [Kj.) The Yukawa couplings of the lightest two families have negligible 
effects on the running of other parameters and will be ignored henceforth. 
We have also ignored in Fig. 2 the small contributions from the r-lepton 
Yukawa coupling. In this approximation, all elements of k evolve identically. 
The parameters are chosen so that m t (the top quark mass) and mn remain 
small enough for perturbation theory to be valid. We take a s (Mz) = 0.12, 
a^ 1 (Mz) = 128, and sm 2 9y/(Mz) = 0.233. The running is significant when 
either m t or mn is large (but especially when m t is large). Large m t and mn 
will cause k to decrease by an order of magnitude or more in running from 
the Planck scale to the weak scale. 

The renormalization of k can alter the predictions of neutrino mass mod- 
els. For example, in SO(10) type theories |y], [12], |l3l , if the neutrino Dirac 
mass matrix is identified with that of the up-quarks and the Majorana mass 
matrix of vr is taken to be identity, the physical neutrino masses will be 
given by the see-saw formula 



ml 

m u ~ — - 
" "r 



M, 



k{M z 



x. 



Y u (m u ) 



(10) 



AMx) 

where m u denotes the physical mass of the up-type quark and M R ~ M x is 
the Majorana mass of z/r. Assuming an intermediate scale Mx ~ 10 10 GeV 
and taking into account the running of Y u (but not the running of k), the 



^ 2 2 2 

authors of Ref. (Tl|] obtain m Ue = 0.05-^-, m Vll = 0.07-p^-, m Vr = 0.18-^. 



Including the running of k reduces these estimates of neutrino masses by 
by an additional factor of 2 to 3 for m^ near 200 GeV. Consequently, 
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the intermediate symmetry SU(3) C x SU(2) L x SU(2) R x U(l) studied in 



Ref. |L2| may become compatible with the Mikheyev-Smirnov-Wolfenstein 



solution of the solar neutrino puzzle. In contrast, large running would 



spoil proposals |15| to relate the just-so solution of the solar neutrino problem 
(where m v ~ 10 -5 eV) with the Planck scale. 

The running of k also affects model predictions of the neutrino mixing 
angles. In the approximation of two neutrino mixing, the mixing angle is 
given by 

tan2fl = 2Kl2 , (11) 

^22 - «11 

in the basis where Y e is diagonal. Eq. (|8|) determines the evolution equation 
for the mixing angle to be 

2 dsm 2 29 . 2 . 2na\/ 2 2^22 + ^11 nn s 

16tt — = sin 20(1 - sin 26){y 2 - y x ) , (12) 

am/i K22 — 

where y\ and 7/2 are the Yukawa couplings of the corresponding charged 

leptons. As expected, the mixing angle does not run when the mixing is 

maximal or zero. However, significant running can occur if the magnitude 

of K22 — «ii is less than or comparable to y 2 T ~ 10~ 4 . (This is the case in 



some pseudo-Dirac [16] neutrino models.) As a result, large mixing at Mx 
can become small mixing at the weak scale, and vice versa. 

Unlike the elements of k, large running of the mixing angle does not 
require a heavy top quark and/or Higgs boson. Examples are shown in Fig. 
3. We have chosen m t (M z ) = 130 GeV and m H (M z ) = 100 GeV. y 2 is 
taken to be the r Yukawa coupling, so Fig. 3 is relevant for e — r and \x — t 
mixing. Curves a) and b) show that the mixing angle can decrease drastically 
from its values at high energies. This can be achieved for neutrino masses 
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of order 1 eV and Am 2 = m\ — m\ of order 10~ 5 to 10~ 6 eV 2 . Curves c) 
and d) illustrate the case in which small mixing at a high scale can become 
large at a lower scale. This requires k 2 2 £ «ii at the high scale, assuming 
v\ > Hi- The running of the mixing angle now depends more sensitively on 
the degeneracy in k. This is because, according to Eq. (H), ^22 decreases 
with decreasing momentum more slowly than ki%. Consequently they will 
cross at some scale, which leads to the resonance curve shown in c). 

3. Two Higgs Doublets. As is customary, we assume that the Higgs 
doublets transform independently under the discrete symmetry 

-> (13) 
for which the most general Higgs potential is 

-A 4 ($ t 1 $ 2 )($ 2 <f 1) - [y ($1$ 2 ) 2 + h.c.) (14) 

The discrete symmetry is to insure that there are no flavor changing neutral 
Higgs couplings in the dimension 4 terms. With this symmetry, each type 
of charged fermions can couple to only one Higgs doublet. We denote the 
doublet which couples to the charged leptons as $1. However, there are now 
four ways to combine two Higgs fields and two neutrino fields, resulting in 
four operators relevant to the evolution of the neutrino mass. These operators 
fall into two classes according to how the product transform under the 
discrete symmetry. Operators where QiQj transform identically will be mixed 
by renormalization. 
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The two operators in which are even are 
C uv = l^l^ll^ie^ + \^l c lll^ l 2^i + h.c. . (15) 

Calculating the class of diagrams shown in Fig. 1, we find the evolution of 

and k^ 22 \ 

16tt 2 ^ = [-3g 2 2 + 2Ai + 2S U ]^ - \[kV\Y}Y 9 ) + Wn) T « (11) ] + 2A^ 22 \ 

w K ( 22 ) 1 
16vr 2 ^- = \-3g 2 2 + 2A 2 + 2S 2 V 22 > + l [^ 22 \Y}Y e ) + (Y}Y e f k^} + 2W 11 ). 

U 111 J^ti 



Here S 17,1 * is the two-Higgs-doublet generalization of Eq. (H). Its precise form 
depends on which charged fermions $ r couples to. For instance, if $i couples 
to the down type quarks and the charged leptons and $2 couples to the up 
type quarks, then 

S 11 = Tr[3Y ] d Y d + Y*Y e ], S 22 = Tr[3Y^Y u \. (17) 

The two operators in which are odd are 

A* = l&Zfy&Mfaji - \^e kl ) + \&Zll^ l 2 e i5 e kl + h.c. (18) 

Here k)^ (and all previous /c's) are symmetric under interchange of the 
generation indices m and n, while ££ 2 -* is antisymmetric. Again calculating 
the class of diagrams shown in Fig. 1, we find the evolution equations 

16vr 2 ^ = [-3 g 2 2 +2X, + 2\, + S}K^-\\^ 2 \Y}Y e ) + {Y}Y^ 

+2[e 2 \YX)-(Yxm 

16vr 2 ^ = l-9g 2 2 +2\ 3 -2\, + S]e 2) + lie 2 \YX) + (YX) T Z il2) ] 

+1^ 12 \yX)-(yX) t ^% (19) 



(16) 



Only the operator contains a neutrino mass term, however the above 
equations show that the two operators are mixed by renormalization and so 
£( 12 ) must also be simultaneously evolved. 

Different neutrino mass models yield different subsets of these four oper- 
ators below the lepton number breaking scale, Mj. For example, the see-saw 
mechanism in some SO(10) models^ may produce and while the 
Zee model || produces and ^ 12 -*. In fact, one can easily construct a 
see-saw type model which leads to any desired subset of the four operators 
by carefully choosing how the different z/r's transform under the discrete 
symmetries. After the operators are produced, they mix under evolution 
according to Eqs. (|16D and (|l^) given above. 

To illustrate the evolution of the neutrino mixing angles in two-Higgs- 
doublet models, we assume that $1 couples to the charged leptons and the 
down type quarks, while $2 couples to the up type quarks. We take the 
neutrino mass matrix to be 

M v = -~v 2 [kWcos 2 /3 + ^ 22 hm 2 p], (20) 

as might occur in a see-saw model. Here Vi/y/2 is the vacuum expectation 
value of $j, tan/3 = v 2 /vi, and v 2 = v \ + v\. The parameter tan/3 runs 
according to 

167r 2 ^^ = tan/3(5 n -5 22 ). (21) 
am/i 

In general, the evolution of the neutrino mixing angles depends on the top 

quark and Higgs boson masses, so large changes during running are possible. 

4 Both even and odd operators can be produced simultaneously if the discrete symmetry 
is broken at Mj. 
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To illustrate a particularly interesting possibility, we have plotted in Fig. 4 
the mixing angle for a two generation system as a function of momentum. 
The evolution equations for the Yukawa and the quartic scalar couplings can 



be found, e.g., in fllQfl . We have chosen the parameters at the weak scale to 
be y t = 1.35, X x = 0.16, A 2 = 1.13, A 3 = -0.0081, A 4 = -0.061, A 5 = -0.011, 
and tan/3 = 1. This corresponds to ra t = 168 GeV and a Higgs boson spec- 
trum of ra H ± = 47 GeV, rap = 26 GeV, = 67 GeV, ran = 186 GeV, 
where H ± , P, h and H stand for the charged, pseudoscalar, and the two 
scalar Higgs bosons. Such a spectrum is phenomenologically acceptable and 
also guarantees that the Higgs potential remains bounded throughout the 
entire momentum range up to Mx, taken here to be 10 14 GeV. Our choice 
of parameters is such that yt and A 2 are near their infrared fixed point values 
while others are not. We then choose the k's to have, in accordance with naive 
expectations, a large hierarchy at Mx'- Kii = 0.05, Kyf 1 = 0.05, iS^ — 1-0 
(in suitable units) and all other elements to be zero. For this parameter 
choice, there occurs at the weak scale a degeneracy between the diagonal 
elements of the neutrino mass matrix. At this "resonance" the mixing is 
maximal. The dashed line in Fig. 4 is the evolution for all the same param- 
eters except = 0.06. Now the resonance occurs at a higher momentum. 
These curves illustrate that small mixing at the high scale can become large 
and even maximal mixing at the weak scale. With suitable choice of the 
parameters, the opposite can also happen, i.e., large mixing at the high scale 
can become small at the weak scale. 

The evolution of the elements of k^ 11 ' and is displayed in Fig. 5. (All 
elements of kS 11 ' evolve identically if the small r Yukawa coupling is ignored, 
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similarly for k^.) Notice that a large variation in (by a factor of 20) 

is now possible even for moderate values of m t . It is clear from Fig. 5 that, 

unlike in the standard model, "resonant mixing" can occur in the two-Higgs- 

doublet models even for nondegenerate neutrinos. The reason is that the 

variation of z^ 22 -* is more prominent than k^ u \ since A2 and y t are large while 

(22) (11) 

Ai is not. As a result, although k 2 2 ^ K ii a t Mx, their values can become 
closer at lower momenta, as can be seen from Fig. 5. The momentum scale at 
which mixing angle resonance occurs does not correspond to the momentum 
at which the two k's cross, since the running of tan/3 also affects the evolution 
of the mixing angle. 

4- Minimal Supersymmetric Standard Model (MSSM). In the MSSM, 
when i?-Parity is assumed to be an exact symmetry, lepton number violation 
can arise only through the dimension 5 superpotential term 

W = ^LTtfHnHveutji • (22) 

Here if 2 is the Y = 1/2 Higgs superfield and L the leptonic doublet superfield. 
Unlike in the non-SUSY two-Higgs-doublet model, the operator in Eq. (22) 
does not mix with any other operators^] The evolution equation for k^ s ' is 
found to be 

+K^Y e +(Y e i Y e ) T «W . (23) 

For a two-family system, the mixing angle evolution is obtained by mul- 
tiplying the right-hand side of Eq. (JTTJ) by (—2). In Fig. 6 we plot the 

5 There are five other dimension 5 operators in the MSSM: (H1H2) 2 , Qu c Le c , 
QQu c d c , QQQL and u c u c d c e c . All these carry zero B — L charge and do not mix with 
the operator of Eq. (22) which has B — L = —2. 
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) tlrv 
ril n / / 



K 



-1g\ - Qg 2 2 + QTr 



Y]Y„ 



running of sin 2 26 l as a function of fi. The parameter choice corresponds to 
Vt — Vb — Vr — 3 at M x = 10 16 GeV, so that tan/3 = m t /m b . The corre- 
sponding top-quark mass is m t = 184 GeV at the weak scale. We choose for 
the solid curve k^] = 1, «i2 =0.1 and «ii = at Mx, which exhibits a large 
hierarchy. sin 2 26 l is seen to increase by about a factor of 2 in running from Mx 
to M z . The dotted curve corresponds to K22 — 1, ^12^ — 

0.035, 0.65 

at Mx- Since y T in the MSSM is larger than in the SM, the running of the 
mixing angle in MSSM is enhanced. Resonant mixing can occur even without 
degenerate neutrinos. 

5. Conclusions. We have calculated the renormalization of the neutrino 
mass operator in the standard model with one and two Higgs doublets, and 
also in the minimal supersymmetric standard model. The renormalization 
group equations are given in Eqs. (|8|), (0), (0) and fl23|). These renormal- 
ization effects should be included in any neutrino mass model. In addition, 
we have solved these equations for certain parameter choices to demonstrate 
some interesting features. If the top quark mass and/or Higgs boson mass are 
large, it is generally true that the neutrino mass operators will have sizeable 
evolution. Their strength can decrease by an order of magnitude in running 
from the Planck scale to the weak scale. Furthermore, running over this 
momentum range can drastically change the neutrino mixing angles; possi- 
bly from small values at high energies to large values at observable scales, 
and vice versa. Large evolution of the neutrino mixing angle requires a near 
degeneracy of the neutrino masses in the standard model, but a degeneracy 
is not required in models with two Higgs doublets. 

Note added. While this work was being completed we received reference 
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[ |j~7|j which discusses the same issues. Our results for the evolution equa- 
tions in the supersymmetric model are in agreement. However, there are 
some differences for the standard model and the two-Higgs-doublet models, 
presumably due to the neglect of our Fig. If in Ref. |T7||. Our analysis is 
more general and reveals the interesting possibility of resonant running of 
the neutrino mixing angles. 
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FIGURE CAPTION 



Fig. 1 One-loop Feynman diagrams contributing to the renormalization of 
k in Landau gauge. 

Fig. 2 Running of k below M x = 10 19 GeV for the standard model with 
one Higgs doublet. The top quark and Higgs boson masses are chosen to be, 
respectively, in GeV: a) 226 and 243.2, b) 220 and 233.8, c) 180 and 183, d) 
130 and 160, e) 130 and 100. 

Fig. 3 Illustration of the running of the neutrino mixing angle in the stan- 
dard model. Parameter choices are: a) Ku(M x ) = 1.0 (in suitable units), 
k 22 {M x ) - K U (M X ) = 10~ 6 , and M x = 7.88 x 10 18 GeV; b) same as curve a) 
except k 22 (M x )-k u (M x ) = lO" 5 ; c) k u (M x ) = 1.0001 and k 22 (M x ) = 1.0; 
d) Kn(Mx-) = 1.0002 and k 22 (M x ) = 1.0 

Fig. 4 Examples of the running of the neutrino mixing angle in two-Higgs- 
doublet models. Parameter choices are given in the text. 

Fig. 5 The running of k* 11 ) and corresponding to the dotted curve in 
Fig. 4. 

Fig. 6 Examples of the evolution of the neutrino mixing angle in the MSSM. 
See text for the choice of parameters. 
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